The low-energy ARPES and heat capacity of Na . 3 CoO 2 : A DMFT study 
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We use the dynamical mean-field theory (DMFT) to calculate the angle resolved photoemission 
spectrum (ARPES) and heat capacity for Nao.3Co02- Both the traditional Hirsch-Fye Quantum 
Monte-Carlo technique and the newly developed continuous time quantum Monte-Carlo technique 
are used to solve the DMFT impurity problem. We show that the e s ' hole pockets on the Fermi 
surface are suppressed as the on-site coulomb repulsion is increased. A quantitative comparison 
with ARPES experiments and bulk heat capacity measurements indicate that the on-site coulomb 
repulsion is large relative to the LDA bandwidth. 



The cobaltatcs have demonstrated a wide variety com- 
plex behavior. The Na rich region of the phase diagram 
displays various degrees of anomalous behavior, such as 
Curie- Weiss behavior near a band insulator [y, charge 
disproportionationQ, and non- Fermi- liquid behavior in 
the resistivity 0| . Alternatively the Na poor region of 
the phase diagram appears to be a Fermi-liquid. The 
magnetic susceptibility displays Pauli behavior, the re- 
sistivity is roughly quadratic at low temperatures [H, and 
the system appears to be homogeneous [2[. Therefore, the 
Na poor region of the phase diagram seems like a natural 
starting point to attempt to explain the ARPES experi- 
ments and heat capacity measurements from a quantita- 
tive standpoint. 

In Naa;Co02, the cubic component of the oxygen crys- 
tal field splits the Co d manifold into a set of 3-fold t2 g 
orbitals and 2-fold e g orbitals, while the trigonal com- 
ponent will further split the t 2g orbitals into a\ g and e' g . 
The nominal valence of Co in this system will be 4 — x, so 
the Fcrmi-cncrgy will fall within the t2 g manifold. The 
LDA band structure displays two degenerate eigenvalues 
and one non-degenerate eigenvalue at the T-point, corre- 
sponding to the e' g and a± g eigenvectors. The splitting 
between the eigenvalues is roughly 1 eV with the e' g levels 
below the Fermi energy and the a\ g above. Despite this 
distinct splitting at the T-point, the on-site orbital en- 
ergies are nearly degenerate. Additionally, the projected 
density-of-states (DOS) clearly show that the a\ g orbital 
character is strongest at the top and bottom of the band 
while the e' g is present through most of the energy range 
of the ti g bands. The Fermi surface consists of a large 
<X\ g pocket around the T-point and six small e' g satellite 
pockets 0. 

Several experimental ARPES studies have been per- 
formed for Nao.3Co02 0, S @, 0]- A general caricature 
of the LDA bands can be seen in the ARPES. The most 
notable difference as compared to LDA is the significant 
narrowing of the bands, and the suppression of the e' g 
pockets below the Fermi energy. Two previous stud- 
ies addressed the effect of correlations on the electronic 
structure for x = 0.3, and they reached completely oppo- 



site conclusions. Zhou et al performed Gutzwillcr calcu- 
lations for a three-band model corresponding to the LDA 
t2 g band structure 0]. Using an infinite on-site coulomb 
repulsion, they show that the quasi-particle bands are sig- 
nificantly narrowed and the e' g hole pockets are pushed 
beneath the Fermi surface. Although the removal of the 
e' g pockets agrees with the ARPES experiments, it is not 
clear if U = oo is an excessive assumption and there- 
fore smaller values of the on-site coulomb repulsion must 
be considered. Ishida et al0] performed DMFT calcula- 
tions for the three-band t2 g states of the cobaltates and 
found that electronic correlations narrow the bands and 
enhance the e' g hole pockets, completely opposite to what 
was found by Zhou et al. 

Singh et al. have proposed that the inclusion of the 
realistic ordering of the Na destroys the e' g hole pockets, 
as demonstrated by LDA calculations for Nan 7Co09[lcj]. 
Given that the pockets will inevitably be destroyed as 
x — > 1 because the chemical potential moves towards the 
top of the t2g bands, Nao.7Co02 is an extremely lib- 
eral test case in which the pockets are barely present in 
the first place. Previous LDA calculations of the realis- 
tic structure which explicitly included the Na, both with 
and without water, had already shown that LDA predicts 
the survival of the pockets for NaiCo02 11, |l2j, even 
when including full structural relaxations There- 
fore, one can safely conclude that the Na potential is not 
the dominant mechanism which destroys the pockets in 
Na . 3 CoO 2 . 

In this study we resolve the issue of the qualitative 
behavior of the pockets. We calculate the ARPES spec- 
trum and the heat capacity at a range of different U 
in order to determine the best agreement with experi- 
ment. In particular, we focus on the presence or absence 
of the pockets, the linear coefficient of the heat capacity 
and the average Fermi velocity measured in ARPES. The 
presence or absence of the pockets will have important 
consequences for the heat capacity given that they arc the 
dominant contribution to the DOS at the Fermi energy, 
and therefore these two issues are intimately connected. 

We perform DMFT calculations for the t2 g bands of the 
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cobaltatcs, represented by the following Hamiltonian: 

+ J2 A ( n a lg ta -n e > girT ) (1) 



where a, (3 are the orbital indices (ie. 



and 



i,j are site indices, a is the spin index, and A is the 
crystal-field splitting between the a\ g and e' g orbitals. We 
use the low-energy hopping parameters t a p and A which 
were fit to the LDA bands by Zhou et al[j|, allow- 
ing for a direct comparison. Ishida et al did not publish 
their hopping parameters, but they appear to be similar 
given the respective bare DOS which were published^. 
We assume the traditional orbital-independent double 
counting [l^ of U(N — |), and therefore the Hartree- 
Fock terms generated by DMFT will be relevant. Be- 
low we will show that A is a key parameter in deter- 
mining the fate of the e' g pockets and therefore the heat 
capacity. Given that LDA is only an approximate tech- 
nique to generate the low energy hopping parameters, 
we will systematically explore the effect of A on the 
results. The LDA value fit by Zhou et al is roughly 
A = —lOmeV, and therefore the on-site orbital energies 
are nearly degenerate. Alternatively, quantum chemistry 
calculations [14J yield a value of 300mey, so one might 
anticipate — lOmeV < A < 300meV. 

DMFT maps the interacting lattice problem onto an 
impurity problem where the non-interacting bath func- 
tion is determined self-consistently [HI]. The effective 
impurity problem is solved using two different impu- 
rity solvers to ensure that the answer is robust. The 
traditional Hirsch-Fye quantum Monte-Carlo (HFQMC) 
method was used (see ref [l3, 15 1 for detailed review), in 
addition to the newly developed continuous time quan- 
tum Monte-Carlo (CTQMC) method CTQMC 
is generally more efficient than the HFQMC and allows 
one to access significantly lower temperatures and larger 
U. 

The e' g pockets are observed in the k-space direction 
corresponding to the real-space direction which connects 
nearest-neighbor Co atoms. ARPES experiments predict 
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the e' g pockets to be below the Fermi energy 
and this is supported by the analysis of the de Haas 
van Alphcn experiments (ljj ]. Given that two previous 
theoretical studies reach opposite conclusions regarding 
the fate of the e' g pockets, we explore this question in 
detail. When determining the Fermi surface from the 
Dyson equation, only the bare Hamiltonian and the self- 
energy at zero frequency are needed [l^|. Within single- 
site DMFT the self-energy is momentum-independent 
and therefore the self-energy at zero frequency acts as 
a renormalization of the on-site e' g and a\ g energy lev- 
els. Therefore, it is useful to reverse-engineer this prob- 
lem and determine what effective e' g and a\ g levels are 



needed to destroy the pockets. This is not a-priori obvi- 
ous given that the neighboring e' g -a\ g hopping is compa- 
rable to the aig-aig hopping. We find that the pockets 
are completely insensitive to perturbations of the a\ g on- 
site energy. The pockets are destroyed when the e' g on- 
site energy is shifted down by roughly 70meV, indepen- 
dent of the value of the perturbation of the a\ g level, and 
therefore the criterion for the destruction of the pockets 
is — fj, < —70meV. This is an important observa- 
tion which indicates that a relatively small perturbation 
of the e s ' on-site energy will destroy the pockets. DMFT 
calculations can now be performed to determine when 
this criterion is satisfied. 
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FIG. 1: a.) e' g orbital occupation as a function of U. The 
dig orbital occupation can be found from the relation n aig = 
5.3 -2n e / HFQMC and CTQMC calculations were performed 

at j3 = 40eV"~ 1 and j3 = 100eV _1 , respectively. Filled points 
indicate that the pockets have been destroyed, b.) The self- 
energy of the e'g orbital at zero frequency minus the chemical 
potential for various values of the crystal field splitting A 
(A = —lOmeV corresponds to the LDA crystal field). 



— fx is found to be a monotonically decreasing func- 
tion of U, which indicates that the pockets are diminished 
with increasing interactions (sec figure [T] b) . The quan- 
tity S ai - /i is a monotonically increasing function of 
U (not pictured). Considering the dynamical portion 
of the self-energy (ie. Yid yn {iw) = £(zo>) — £(oo)), we 

find that £^ n (0) > 2^(0), but this is countered by 
the fact that £ ai9 (oo) > £ e s(oo) and the net result is 
that S Q1 «(0) > E e s(0). Therefore, the destruction of the 
pockets is due entirely to the static Hartree contribution 
to the self-energy (ie. £(oo) ). The dynamical contri- 
bution of the self-energy actually enhances the pockets, 
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but ultimately the static part of the self-energy domi- 
nates and the net effect is that increasing interactions 
diminishes the pockets. The preceding analysis is true 
for both the CTQMC and the HFQMC impurity solvers. 
The two solvers agree completely on a qualitative level, 
and are quantitatively similar for small U while more sig- 
nificant differences arise for larger U due to discretization 
errors within the HF method. We also plot the orbital 
occupations as a function of U (see figure QJi), which in- 
dicates an enhancement of orbital polarization as U is 
increased, consistent with the destruction of the pock- 
ets. This conclusion is found for both the CTQMC and 
HFQMC solvers, and this holds even when an on-site 
exchange coupling (ie. —Jn aa n^ a ) is included in the 
HFQMC method. Exchange had little effect on the re- 
sults, in agreement with Ishida et al [§]. 

Having established the qualitative behavior of the 
pockets, we continue with more quantitative analysis. 
Given that only a 70 meV downward shift of the e g ' 
on-site energy is required to destroy the pockets, we be- 
lieve that it is useful to probe the behavior of the self- 
energy for other values of A as rationalized above. The 
CTQMC calculations were repeated for A = 50meV and 
A = lOOmeV. As anticipated, starting with a crystal 
field splitting which diminishes the pockets acts cooper- 
atively with interactions and causes the system to polar- 
ize more and the pockets to be destroyed for a smaller 
U (see figure Q] a and b). For A = 50meV the pockets 
are destroyed for U > 2.0eV, while for A = lOOmeV the 
pockets are destroyed for U > l.OeV. Therefore, we con- 
clude that the value of U required to destroy the pockets 
depends strongly on the value of A. 

Our findings are in qualitative agreement with Zhou 
et al[|| and in qualitative disagreement with Ishida et 
aid. It is not clear why the results of Ishida et al are 
different, but it is likely the result of differences in the 
bare Hamiltonian. Above we showed that the static cor- 
relations diminish the pockets while the dynamic corre- 
lations enhance the pockets, and the balance of these 
two effects will be influenced by the bare Hamiltonian. 
Regardless, our following analysis of the Fermi velocity 
and heat capacity will demonstrate that the destruction 
of the pockets is the only viable possibility to achieve 
agreement with experiment. 

Having understood the behavior of the pockets, we 
now compute the linear coefficient of the heat capacity, 
7. We use Fermi-liquid theory to calculate 7 from the 
DOS at the Fermi energy and the quasiparticle weight 
Z which is calculated in our QMC calculations. Math- 

2-nk 2 R v-^ Pa(0) 1 

r v where a corre- 



cmatically, we have 7 
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sponds to the orbital index and p is the local spectral 
function [l3j|. The experimentally measured heat capaci- 
ties for Nan^CoO^ are found to be in the range of 12-16 

mj 
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FIG. 2: The heat capacity, density-of-states at the Fermi 
energy, the quasiparticle weight Z, and the average of the 
absolute value of the Fermi velocity for the T-M and T-K 
directions are shown in panels a, b, c, and d, respectively. 
Filled points indicate that the pockets have been destroyed. 
The shaded region in panel a corresponds to the range of 
experimentally measured values of 7. The blue dotted line 
in panel d corresponds to the experimentally measured Fermi 
velocity, while the grey shading corresponds to the error bar 
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hoppings, which is already within the bounds of experi- 
mental measurements. This may mislead one to believe 
that correlations are negligible, but a more careful exam- 
ination shows otherwise. The DOS at the Fermi energy 
initially decreases weakly as U increases and eventually 
drops in a discontinuous fashion, which signifies the de- 
struction of the pockets (see figure [2]). The quasiparticle 
weight also decreases as U increases. Given that the lin- 
ear coefficient of the heat capacity is proportional to the 



ratio 



p(e/: 



22, 23. 24, 2i 



mol-Co-K 2 

that the heat capacity is 14.27 



We begin by noting 
for the LDA 



mol-Co-K 2 



, the overall effect is not apriori obvious. The 
heat capacity initially increases as U increases, then dis- 
continuously drops when the pockets are destroyed, and 
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eventually plateaus for large U. Increasing A causes the 
drop in the heat capacity to occur at smaller values of 
U and a overall lower value for the heat capacity. In 
order to achieve agreement with experimental measure- 
ments of the heat capacity, one needs a relatively large 
U = 6eV when using the LDA A, and even larger U is 
needed for larger A. A key point is that if the pockets are 
retained, the 7 becomes excessively large as U increases. 
Given that the ARPES indicates that the bandwidth is 
nearly halved as compared to LDA, an appreciable U 
must be present to narrow the LDA bands and if the 
pockets were still present the heat capacity would be ex- 
cessive as compared to experiment. It is reasonable to 
expect that the heat capacity should be under-predicted 
when only considering the Hubbard model. There will 
likely be electron-phonon coupling to the local breathing 
mode of the octahedron, or perhaps other modes, which 
will induce a narrowing of the bands and therefore an 
enhancement of the heat capacity. 

The experimentally measured Fermi velocity may also 
be calculated as a function of U (see figure [2]). Increasing 
the U decreases the quasiparticle weight Z and therefore 
decreases the Fermi velocity. In order to achieve veloci- 
ties comparable with experiment, one needs a relatively 
large U > 3eV. This is another piece of evidence, inde- 
pendent of the heat capacity measurements, which indi- 
cates that the U must be relatively large. Once again, if 
the U is large then the pockets must be absent in order 
to get acceptable agreement with the heat capacity. 

Both impurity solvers used in this study work on the 
imaginary axis, and therefore one must perform an ana- 
lytic continuation to access real frequency quantities like 
the ARPES spectrum. Various approaches exist to per- 
form the analytic continuation, but all arc approximate. 
We expand the self-energy to first-order, which allows an 
exact analytic continuation, and use the resulting self- 
energy to construct the low energy ARPES spectrum (see 
figure[3|). The first case corresponds to the A = —lOmeV 
and U = 6eV, the minimum to destroy the pockets for 
this case (see figure As shown, the pockets are 

too close to the Fermi energy, and the bands are exces- 
sively narrow as compared to experiment. Increasing to 
U = lOeV will further sink the pockets, but it will also 
to continue to narrow the bandwidth. Alternatively, one 
may decrease the U slightly if the A is increased as the 
pockets are suppressed at a much lower U . Therefore, 
we examine the ARPES for A = 50meV" and U = 4eV, 
and A = lOOmeV and U = 4eV (see figure [3Jd and[3fc). 
Decreasing U and increasing A increases the bandwidth 
and pulls the pockets down further, putting the result 
in better agreement with experiment. This analysis sug- 
gests that optimum U should be chosen from the lower 
end of the range of values deduced from the analysis of 
the heat capacity and the velocity, and that A should be 
slightly larger than the one deduced from LDA. 

In conclusion, we have examined the issue of the e' g 
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FIG. 3: The ARPES for a.) U = 6.0 A = -lOmeV b.) U = 
4.0 A = 50meV c.) U = 4.0 A = WOmeV. All calculations 
performed at /3 = 100eV -1 . 



pockets, the value of the linear coefficient of the heat ca- 
pacity, and the Fermi velocity. We have demonstrated 
that increasing interactions destabilize the e' g pockets 
and pushes them beneath the Fermi energy. This is in 
agreement with previous calculations of Zhou et alQ and 
in disagreement with previous calculations of Ishida et 
alQ. Reasonable agreement can be achieved with both 
bulk heat capacity measurements and the Fermi veloc- 
ity measured by ARPES when using an on-site coulomb 
repulsion which is several times the LDA bandwidth (ie. 
U > 4eV). 
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